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The ground state phase diagram of the ahernating spin-1/2 chains with anisotropic ferro- 
magnetic coupling under the influence of a symmetry breaking transverse magnetic field is 
studied. We have used the exact diagonalization technique. In the limit where the antiferro- 
magnetic coupling is dominant, we have identified two Ising-type quantum phase transitions. 
We have calculated two critical fields h^ and hc2 , corresponding to the transition between 
different magnetic phases of the system. It is found that the intermediate state (/ici < h < hc2) 
is gapful, describing the stripe-antiferromagnetic phase. 

KEYWORDS: Alternating Heisenberg Chains, Transverse Magnetic Field 



1. Introduction 

The effect induced by external magnetic fields in the 
low-dimensional magnets has attracted much interest re- 
cently from experimental and theoretical points of view. 
In particular, critical properties of the alternating spin- 
1/2 chains in a magnetic field have been a field of in- 
tense studies. This seems pertinent in the face of great 
progress made within the last years in fabrication of such 
AF-F compounds. Since, the AF-F chains have a gap in 
the spin excitation spectrum, they reveal extremely rich 
quantum behavior in the presence of the magnetic field. 

A typical example of the ferromagnetic-dominant AF- 
F bond alternating chain is (CH3)2CHNH3CuCl3 (iso- 
propylammonium coppertrichloride: IFACuCla).^'^' The 
energy gap in the absence of the external magnetic field 
is estimated from the susceptibility to be 17 — ISK.-'^' 
This value is also confirmed by the analysis of the spe- 
cific heat.^^ From the viewpoint of the crystal structure, 
the origin of the spin gap was expected to be the spin- 
1/ 2 AF-F alternating chain along the c-axis.^^ However, 
quite recently, it was suggested that this system should 
be characterized as a spin ladder along the a-axis with 
the strongly coupled ferromagnetic rungs, namely the an- 
tiferromagnetic chain with effective S = 1, and the ex- 
citation gap A was estimated as 13. 6K by means of the 
neutron inelastic scattering experiments.^' 

More recently, it was reported new inelastic neutron 
scattering results for (CH3)2NH2CuCl3 (Dimethylammo- 
nium copper II chloride, also known as: DMACUCI3 or 
MCCL). The linked-clustcr calculations and the bulk 
measurements show that DMACUCI3 is a realization of 
the spin-1/2 alternating AF-F chain by nearly the same 
strength of the antiferromagnetic and ferromagnetic cou- 
plings.**' 

There are other examples for the AF-F alternating 
spin-1/2 chain compounds such as: [Cu(TIM)]CuCl4 
(TIM=2,3,9,10-tetramethyl-l,3,8,10-tetraenecyclo- 
1,4,8,11-tetraazatetradecane) and (4 — BzpipdH)CuCl3 
(4-BzpipdH=4-benzylpiperidinium) ,''^' ^' and 

CuNb20&.'^^ 



Theoretically, the AF-F alternating chain is expected 
to have relation to the Haldane-gap systems,^-' since they 
are regarded as the spin-1 antiferromagnetic chain in the 
large ferromagnetic coupling limit. The energy gap exist 
between the ground state and the first excited state. The 
spin correlation function of the ground state decays ex- 
ponentially as well as the spin-1 AF chain. The ground 
state properties^"'^^^ and low-lying excitations'^^^ of this 
model were well investigated by numerical tools and vari- 
ational schemes. In particular, the string order parame- 
ter originally defined for the spin-1 Heisenberg chains^^' 
was generalized to this system. The ground state has 
the long-range string order, which is characteristic of the 
Haldane-gap phase. Hida has shown that the Haldane 
phase of the AF-F alternating chain is stable against any 
strength of the randomness.*^' 

The ground state phase diagram of the AF-F alternat- 
ing chain in a magnetic field is studied by the numerical 
diagonalization and the finite-size scaling based on the 
conformal field theory.*^) It is shown that the magnetic 
state is gapless and described by the Luttinger liquid 
phase. It is also found that the magnetic state is char- 
acterized by the algebraic decay of the spin correlation 
functions. Recently, Yamamoto et. al described the mag- 
netic properties of the model in a magnetic field in terms 
of the spinless fermions and the spin waves. *^' They em- 
ployed the Jordan- Wigner transformation and treated 
the fermionic Hamiltonian within the Hartree-Fock ap- 
proximation. They have also implemented the modified 
spin wave theory to calculate the thermodynamic func- 
tions as the specific heat and the magnetic susceptibility. 

The effect of an uniform transverse magnetic field 
on the ground state phase diagram of a spin-1/2 AF- 
F chain with anisotropic ferromagnetic coupling is much 
less studied. Partly, the reason is that AF-F alternating 
chains with anisotropic ferromagnetic coupling are not 
still fabricated. However, from the theoretical point of 
view these systems are extremely interesting, since they 
open a new wide polygon for the study of complicated 
quantum behavior, unexpected in the more conventional 
spin systems. The Hamiltonian of the model under con- 
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Fig. 1. Sketch of the Spin-1/2 AF-F alternating Chain considered 
in this paper. 



sideration on a periodic chain of N sites is given by 

N/2 
H = JaF 2_^[S2j^lS2j + S2j^iS2j + S2j^iS2j] 



+ 



N/2 
Jp / ,['-'2 j ■-'27 + 1 

N 



qy qV 
'^2j^2j+l 



^^2j^2j + l\ 



(1) 



where S^'^'^ are spin-1/2 operators on the j-th site. Jp 
and Jaf denote the ferromagnetic and antiferromagnetic 
couphngs respectively. The limiting case of isotropic fer- 
romagnetic couphng corresponds to A = 1 and h is the 
transverse magnetic field. The model (1) in the case of 
Jp = —Jaf and A = 1 corresponds to the isotropic 
Heisenberg chain in an external magnetic field. The 
ground state phase diagram of this model is known. ^°^ 

In this paper, we present our numerical results ob- 
tained in the low-energy states of the AF-F alternating 
spin-1/2 chain with anisotropic ferromagnetic coupling 
in a transverse magnetic field h. Assuming that the an- 
tiferromagnetic coupling is dominant, we study the ef- 
fect of a uniform transverse magnetic field on the ground 
state phase diagram of the model. In particular, we ap- 
ply the modified Lanczos method to diagonalize numer- 
ically finite chains. Using the exact diagonalization re- 
sults, we calculate the spin gap, the magnetization, the 
string order parameter, and various spin-structure fac- 
tors as a function of the uniform transverse magnetic 
field. Based on the exact diagonalization results we ob- 
tain the magnetic phase diagram of the model show- 
ing Haldane, stripe-antiferromagnetic and ferromagnetic 
phases. We denote by "ferromagnetic phase" the phase 
with the magnetization parallel to the external field as 
only the nonvanishing order parameter. We show that 
the Haldane phase remains stable even in the presence 
of the uniform transverse field less than a critical field 

The outline of the paper is as follows: In section 11 we 
discuss the model in the strong antiferromagnetic cou- 
pling limit and derive the effective spin chain Hamilto- 
nian to outline the symmetry aspects of the problem. In 
section III we present results of the exact diagonalization 
calculations using the modified Lanczos method. Finally 
we conclude and summarize our results in section IV. 



2. Large Antiferromagnetic Coupling Limit 

In this section we will discuss the model (1) in the 
limiting case of the strong antiferromagnetic coupling 
Jaf ^ Jf- We will show that in this limit the model 
can be regarded as a fully anisotropic XYZ chain, which 
allows us to outline the symmetry aspects of the problem 
under consideration.^^) The schematic picture of the AF- 
F alternating chain is plotted in Fig. 1. At Jaf S> Jf, 
the system behaves as the nearly independent block of 
pairs. ^^) Indeed an individual block may be in a singlet 
(|5 >= l/\/2[| ti> -| iT>]) or a triplet state (jTi >= 
I TT>, \To >= 1/V2[| Ti> +1 iT>], \T-i >= I ii>) with 
the corresponding energies given by 

3 Jaf 



Ei-i — 



Jaf _^ , \ J-, Jaf 



E, = — 



At ft. < Jaf, one component of the triplet becomes closer 
to the singlet ground state such that for a strong enough 
magnetic field we have a situation when the singlet and 
5*^ = 1 component of the triplet create a new effective 
spin T — 1/2 system. One can easily project the original 
Hamiltonian (1) on the new singlet-triplet subspace 



1 



i^)^i.) = -^[in)-uT)], 



I t) ^ |Ti) = I tt) • 



When expressed in terms of the effective spin operators 
Tj up to the accuracy of irrelevant constant, the orig- 
inal Hamiltonian (1) becomes the Hamiltonian of the 
spin-1/2 chain fully anisotropic XYZ chain in an effective 
magnetic field^^' ^^^ 



H. 



eff 



N/2 



^1 ^1+1 



-f-]Vi + ATfrlVJ 



JV/2 



,eff 






(2) 



where h'^f^ = h — Jaf + Jp/^- Note that in deriving 
(2), we have used the rotation in the effective spin space 
which interchanges the x and z axes. 

At A = 1, the effective problem reduces to the the- 
ory of the XXZ chain with a fixed antiferromagneic 
anisotropy of 1/2 in a magnetic field. The gapped phase 
at h'^^f < h%l^ for the AF-F alternating chain corre- 
sponds to the negatively saturated magnetization phase 
for the effective spin chain, whereas the massless phase 
for the AF-F alternating chain corresponds to the finite 
magnetization phase of the effective spin-1/2 chain. The 
critical field h'^JJ where the AF-F alternating chain is 
totally magnetized, corresponds to the fully magnetized 
phase of the effective spin-1/2 chain. 

From the exact ground state phase diagram of the 
anisotropic XXZ chain in a magnetic field, ^''^ it is easy 
to check that ftc{c2 ~ TJf/4. Therefore we can see 
that the isotropic AF-F alternating chain in a magnetic 
field shows a transition from the Haldane phase to the 
Luttinger-Liquid phase at h^ = Jaf — Jf/'^ and a tran- 
sition from the Luttinger-Liquid phase into the fully po- 
larized phase at h%^* = Jaf- 

Away from the isotropic point A = 1 the effective 
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Fig. 2. Difference between the energy of the two lowest levels and 
the ground state energy as a function of the transverse magnetic 
field, h for J ap = 2Jf and A = 0.5 including different chain 
lengths A'' = 12, 16, 20. The lowest three lines for E\ — Eq arc 
indistinguishable on this scale. 



Hamiltonian (2) describes the fully anisotropic ferromag- 
netic XYZ chain in a magnetic field that is directed per- 
pendicular to the easy axes. For the particular value of 
magnetic field h'^^^ = 0, the effective XYZ chain is long 
range ordered in "y"— direction. This is corresponding 
to the original AF-F alternating chain being ordered in 
the direction perpendicular to the applied magnetic field 
with opposite staggered magnetization on the blocks. We 
have defined this new phase as a stripe-antifcrromagnetic 
phase (Fig. 6). The stripe-antiferromagnetic phase ex- 
tends over the whole region of the transverse magnetic 
field ft-d < h < hc2- For other values of the h'^^^ it is 
clear that the stripe-antiferromagnetic phase order will 
be replaced either by the Haldane phase {h < hc^ corre- 
spond to negative effective fields) or the phase with only 
one order parameter - magnetization along the applied 
field {h > hc2 correspond to positive effective fields). 

3. Numerical Results 

In this section, for exploring the nature of the spec- 
trum and the phase transition, we use the modified 
Lanczos method to diagonalize numerically finite chains 
{N — 12, 16, 20, 24). The energies of the few lowest cigcn- 
states were obtained for the chains with periodic bound- 
ary conditions. To find the effect of a transverse magnetic 
field on the ground state phase diagram of the system, we 
start our consideration by the anisotropic case, A ^ 1. 
First, we have computed the three lowest energy eigenval- 
ues of A^ = 12, 16, 20 chain with Jaf = "^Jp and different 
values of the anisotropy parameter A — 0.25, 0.5, 0.75. As 
an example, in Fig. 2 we have plotted results of this calcu- 
lations for A = 0.5. It can be seen that the difference be- 
tween the two lowest states, decreases by increasing the 
transverse magnetic field. These are independent of the 
chain length. Considering the limiting case oi Jp = the 
system consists of the N/2 independent block of pairs. 
Indeed the energy of the ground state for the A^/2 blocks 



Fig. 3. The transverse magnetization M^ as a function of applied 
field h for A'^ = 20 chain with Jaf = 2 Jf and for different values 
of the anisotropy parameter A = 0.25, 0.50, 0.75. 



is equal to A^/2 x Eg, and in the presence of a mag- 
netic field the first excited state has the energy equal to: 
(A^— 1)/2 X Es+Ei. Thus the deference between two low- 
est energy of the system will be independent of the size 
of the system. However for the case of the Jap ^ Jp the 
other terms in the Hamiltonian behave as a small per- 
turbation, which causes that the three lines for Ei — Eq 
to be indistinguishable on this scale. 

In the region of the magnetic fields /ici < h < hc2 
the two lowest states form a twofold degenerate ground 
state in the thermodynamic limit. Thus, the excitation 
gap in the system is the difference between the second 
excited state and the ground state. In the case of the 
zero magnetic field {h = 0) the spectrum of the model is 
gapped. The spectrum remains gapfuU except at the two 



critical fields: 



he, = 1.64 ±0.01 and h^^ = 1.88 ± 0.01. 



The spin gap, which appears at ft, > h^ , first increases 
vs external field, but then starts to decrease, and finally 
vanishes at ftca . At /i > hc2 the gap opens again and, for 
a sufficiently large field becomes proportional to h. These 
results are in good agreement with the results obtained 
in the studies of the fully anisotropic antifcrromagnetic 
XYZ chain in a magnetic field. ^^^ 

On the other hand, we have also checked our numer- 
ical tools in the isotropic case, A = 1. In this case, we 
found that the gap decreases linearly and vanishes in the 
critical magnetic fields h^ = 1.57 ± 0.01 and hc2 =■ 2.0, 
for the constant coupling Jap = 2Jp. Both values of 
the critical fields are obtained from studying of the fi- 
nite chains are very close to the previous section results. 
Also, our numerical results showed that the spectrum re- 
mains gapless for he, < h < hc2 , in good agreement with 
numerical results reported by Sakai.^^^ 

To recognize the different phases induced by the trans- 
verse magnetic field in the ground state phase diagram, 
we have implemented the modified Lanczos algorithm 
of finite size chains (A^ = 12, 16, 20, 24) to calculate the 
order parameters and the various spin correlation tunc- 
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Fig. 4. The string correlation function Ostr{l, ^) for the I = N/2 
as a function of the transverse magnetic field, h for Jaf = 2Jp 
and A = 0.5 including diff'erent chain lengths A'^ = 12, 16, 20, 24. 
The inset has shown the string order parameter, Ostr{N/2, N), 
as a function of the 1/N for a value of magnetic field h = 1.7 > 



Fig. 5. Fixed points /ici (N, N + 4) and h^^ {N, TV + 4) of the phe- 
nomenological renormalization group equation for the energy 
gap, plotted versus 1/{N + 2), for Jaf = 2Jf and A = 0.5 
including diff'erent chain lengths N = 8, 12, 16, 20, 24. The ex- 
trapolated values arc: feci = 1.64 it 0.01 and hc2 = 1-88 ± 0.01. 



tions. 

The first insight into the nature of the different phases 
can be obtained by studying the magnetization along the 
field axis 



^^ = ^E(^^) 



(3) 



where the notation (...) represent the expectation value 
at the lowest energy state. In Fig. 3 we have plotted 
the magnetization along the applied transverse magnetic 
field, M^, vs h for the chain of length N — 20. For arriv- 
ing at this plot we considered Jaf = 2Jf for the different 
values of the anisotropy parameter A = 0.25,0.5,0.75. 
Due to the profound effect of the quantum fluctuations 
the transverse magnetization remains small but finite for 
< h < hci and reaches zero at h = 0. This is in agree- 
ment with the results of magnetization obtained in the 
case of fully anisotropic XYZ chain. ^^^ Increasing the 
magnetic field, the magnetization increases for h > hc^ 
very fast. This behavior is in agreement with expecta- 
tions, based on general statement that in the gapped 
phase, the magnetization along the applied field appears 
only at a finite critical value of the magnetic field equal 
to the spin gap. However, in finite systems we do not 
observe a sharp transition close to the saturation value, 
which happens at /i > hc2 ■ The values of the critical fields 
ft-ci and /ic2 depend on the anisotropy parameter A. By 
increasing A, the critical fields take larger values. Also, 
our numerical results show that the magnetization along 
the directions perpendicular to the applied field (M^ and 
M^) remains zero. 

We employ the phenomenological renormalization 
group (PRG) method^^) to determine these critical fields 
(ft-ci, and hc2)- In the disordered phase [hc-^ > h or 
hc2 < h), the gap {Ei — Eq) tend to a finite value so that 
N{Ei—Eo) increases with N. In the Ising phase (Ei—Eq) 



decreases exponentially with N so that N{Ei — Eq) also 
decreases with N. The phenomenological renormaliza- 
tion group equation is: 



{N + AjgiN + 4, h') ^ Ng{N, h). 



(4) 



Where g{N,h) = Ei{N,h) - Eo{N,h) is the gap value 
for chain length N in the magnetic filed h. At the crit- 
ical point, N{Ei — Eq) should be size independent for 
large enough systems in which the contribution from ir- 
relevant operators is negligible. Due to this situation, 
we can accurately determine the critical points by the 
PRG method. We can define hc{N,N + 4) as the iV- 
dependent fixed point of Eq. 4 and it is extrapolated to 
the thermodynamic limit, in order to estimate he. At 
the critical point h = h' = he, therefore the curves of 
N{Ei — Eq) vs. h for sizes N and A^ -I- 4 cross at cer- 
tain values /ici (A^, A" -I- 4) and hc2 {N, A^ + 4) (finite size 
critical points) . The thermodynamic critical points {hc-^ , 
and /1C2) can be obtained by appropriately extrapolat- 
ing hc,{N,N + 4) or hc2{N,N + A) to N ^ 00. Fig- 
ure (5) represents the extrapolation procedure of the 
transition points, for Jaf = '2Jf and A = 0.5 includ- 
ing different chain lengths A^ = 8, 12, 16, 20, 24. The val- 
ues of hci,2AN,N + 4) for four pairs of system sizes 
{N,N + 4) = (8, 12), (12, 16), (16, 20), and (20,24) are 
represented by square (circle) in this figure. The extrap- 
olated values are: ft-d = 1.64±0.01 and hc2 — 1.88±0.01. 
The system size dependence of the critical points is al- 
most negligible as shown in Fig. 5. 

As we mentioned before, in the absence of the uniform 
transverse field, the spectrum is gapful. The ground state 
of the model is the Haldane phase with the long-range 
string order. Thus, the Haldane phase can be recognized 
from studying the string correlation function. The string 
correlation function in the chain of length A^ defined only 
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Fig. 6. Classically, the long-range order canted in the direction 
of applied field h, which is expected in the ground state of the 
spin-1/2 AF-F alternating chain in a transverse field. 



for odd I as^^^ 



I 2J+1+1 \ 

exp{i7T Y^ S^} 

\ k=2j + l I 



(5) 



In particular we have calculated the string correlation 
function Ostr{l,N) for the I = N/2. In Fig. 4 we have 
plotted Ostr(N/2,N) as a function of h for the chain 
with Jaf — '^Jf, a = 0.5 and for different values of 
the chain length TV = 12, 16, 20, 24. As it can be seen 
from this figure, at /i < /ici, Ostr{N/2,N) is close to its 
maximum value 1.0, therefore the chain system is in the 
Haldane phase. The Haldane phase remains stable even 
in the presence of the transverse field less than hc^ . 

In the inset of the Fig. 4 we have also plotted the 
string order parameter, Ostr{N/2,N), as a function of 
the 1/A^ for a value of magnetic field h = 1.7 > /id. 
It is clear that by increasing the size of the system, the 
Ostr{N/2, N), converges to the very small values close to 
zero. Which shows that there is not the string ordering 
at larger transverse magnetic fields h > hc^ . 

To get additional data about the character of the spin 
ordering in the intermediate gapped phase, we intro- 
duce a new order parameter. Classically, the effect of the 
transverse magnetic field h is interesting. Energetically, 
the ground state of the system has a long-range order 
canted in the direction of the applied transverse magnetic 
field, as illustrated in Fig. 6. The ordering of this phase 
is a kind of the stripe-antiferromagnetic phase. There- 
fore the order parameter of the stripe-antiferromagnetic 
phase has been defined by 



M-- 



^(^1.- 



Sl){' 



1) 



(6) 



For any value of the transverse field h in the interme- 
diate phase, the Lanczos results lead to M^^ = 0, since 
the ground state is two-fold degenerate and in a finite 
system no symmetry breaking happens . However the 
spin correlation function diverges in the ordered phase 
as N — y oo. We have computed the correlation function 
of the stripe-antiferromagnetic order parameter given by 



lN/2 



X" 



~ \ Z^('^2j-1 ^2])\^2j-l+2n ^2]+27i)i ^) 



In Fig. 7 we have plotted x^^ a-s a function of h for 
the chain with Jaf = '^Jf, A = 0.5 and for differ- 
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Fig. 7. The correlation function of the stripe-antiferromagnetic 
order parameter x^^ ^^ ^ function of the transverse magnetic 
field, h for Jaf = 2Jj? and A = 0.5 including different chain 
lengths A'^ = 12, 16, 20, 24. In the inset the correlation function 
of the stripe-antiferromagnetic order x^^ is plotted versus N for 
a value of the transverse field hci < h = l.S < he., ■ 



ent values of the chain length N = 12,16,20,24. As 
it is clearly seen from this figure, there is no long rang 
stripe-antiferromagnetic order along the "y" direction at 
h<hc^ = 1.65±0.05 and h > hc^ ^ 1.90±0.05. However, 
in the intermediate region hc^ < h < hc^ , spins of each 
block show a profound antiferromagnetic order in the 
"y" direction. Thus, between the critical fields /ici and 
hc2 , the ground state of the anisotropic case is Ising-type 
symmetry broken phase (stripe-antiferromagnetic phase) 
unlike the Luttinger liquid state for the isotropic case."'^^-' 
In the inset of Fig. 7 we have also shown the correlation 
function of the stripe-antiferromagnetic order x^^ versus 



TV for a value of transverse field hc^ < /i = 1.8 



< hc2- It 



can be seen that the function x^^ shows the linear de- 
pendence on the chain length N, corresponding to the 
ordered stripe-antiferromagnetic phase in the intermedi- 
ate region hci < h < hc^ ■ 

Thus, our numerical results show that the ground state 
phase diagram of the antiferromagnetic dominant spin- 
1/2 AF-F chain in a transverse field contains, besides the 
gapped Haldane and ferromagnetic phases, the stripe- 
antiferromagnetic phase. Each phase is characterized by 
its own type of long-range order: the ferromagnetic order 
along the transverse magnetic field axis in the ferromag- 
netic phase; the string order along the "z" axis in the 
Haldane phase; and the stripe-antiferromagnetic order 
along the "y" axis in the stripe-antiferromagnetic phase. 

4. Conclusions 

In this paper, we have investigated the ground 
state phase diagram of the antiferromagnetic dominant 
(jiJAF > Jf) spin-1/2 AF-F chain with anisotropic ferro- 
magnetic coupling in a transverse magnetic field h. We 
have implemented the modified Lanczos method to diag- 
onalize numerically finite chains. Using the exact diago- 
nalization results, we have calculated the various order 
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parameters and the spin-structure factors as a function 
of the transverse magnetic field h. We have found that 
a gapped phase exists in the intermediate values of the 
transverse field (/ici < h < hc^)- Then, we have iden- 
tified its ordering as an interesting kind of the stripe- 
antiferromagnetic phase. Two quantum phase transitions 
in the ground state of the system with increasing trans- 
verse magnetic field have been identified. The first tran- 
sition corresponds to the transition from the gapped Hal- 
dane phase to the gapped stripe-antiferromagnetic phase. 
The other one is the transition from the gapped stripe- 
antiferromagnetic phase into the fully polarized ferro- 
magnetic phase. 

On the other hand, in the limit that the antiferromag- 
netic coupling is dominant {Jaf S> Ji^) we have mapped 
the model (1), onto an effective XYZ Heisenberg chain in 
an external effective field (h^''). This mapping allowed 
us to relate the critical fields ft-d and hc2 to the coupling 
constants in the isotropic case A = 1. 
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